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1 Introduction 



The isolation and investigation of the polarized gluon parton distribution is one of the 
most interesting topics for deep inelastic scattering. The three main approaches are 
the analysis of dijet events, on which we will concentrate here, the Q 2 dependence of 
structure functions and the production of charmed particles [0. Up-to-date parton 
distributions based on all experimental input available do not yet determine the form 
of the polarized gluon distribution with any degree of reliability ||. This underscores 
the great significance of a polarized HERA, where all three approaches described above 
would be feasible. 

In the analysis of dijet events two different processes have to be taken into account, 
the photon-gluon fusion (PGF) and the gluon bremsstrahlung (GB), see Fig. |I|. Only 
the photon-gluon fusion allows to analyze the gluon parton distribution, while the gluon 
bremsstrahlung is a background process. In the following we will investigate to which 
extend we can get information about the PGF events from the polarized dijet cross 
section. For this purpose we discuss in the first section, after some remarks to the kine- 
matic cuts used, the influence of quark and remnant masses on the dijet cross section 
in different MC approaches. In the second section we analyze the spin asymmetry of 
dijet events as a function of several kinematic variables. We present a method of how to 
find suitable cuts in these kinematic variables so that the significance of the asymmetry 
signal becomes largest. This is a continuation of the discussion started in ||, where 
the discussion was concentrated on the kinematic variable x g , the fractional momentum 
carried by the gluon. 

The discussion is done on the level of parton jets because on this level the QCD im- 
proved parton model provides a well defined theory, where comparison of different MC 
approaches is possible in a quantitative way. Moreover, this will be the bases for more 
parameter-dependent additions like parton showering, to describe higher-order effects, 
and fragmentation. We should stress, however, that the advances in jet reconstruction 
justify a discussion on this level. For the simulations we use the following two MC 
generators Pepsi and Mepjet. Pepsi is a polarized add-on package to the unpolar- 
ized Lepto-6.5 code. The general principles are the same as in the earlier version 
Pepsi 1 |J. The code contains massless partonic cross sections, where the effects of 
heavy quark masses and remnant target masses are taken into account by constraints 
reflecting energy-momentum conservation of massive particles. 

As to the Mepjet code we will here use two versions. The first one neglects quark 
and target remnant masses, even charm and bottom quarks are treated as massless. 
The second one contains the quark masses in the photon-gluon fusion cross section, but 
neglects the influence of the target remnant mass. We use this fact to demonstrate the 
quark and remnant mass threshold effects in our simulations by comparing the PEPSI 
with the Mepjet results. We show that those effects are small for QCD bremsstrahlung 
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and large for photon-gluon fusion. 



2 Mass effects and kinematics 

The calculations are done for HERA energies, i.e., 820 GeV protons are scattered off 
27.5 GeV electrons. The cuts used are compatible with the present detectors at HERA. 
The option of polarized protons is discussed as a possible HERA upgrade ||. We use 
the following maximal kinematic cuts: 

1(T 5 < x < 1 

0.3 < y < 0.8 (1) 
5 GeV 2 < Q 2 < 100 GeV 2 . 

Q 2 = —q 2 is the virtuality of the photon, x = Q 2 / (2P ■ q) the usual Bjorken x, y = 
(P ■ q)/(P ■ k), and k the momentum of the incoming lepton. The cut in Q 2 is guided 
by the idea that we look for such a kinematic range, where the fusion process is by 
far dominating the bremsstrahlung process (see also [Q). This is, of course an indirect 
effect because the low Q 2 implies small x, which in principle accounts for the large gluon 
contribution. The cut y > 0.3 does not originate from detector requirements, but is due 
to the fact that the spin-dependent part of the cross section, which is proportional to 
y(l — y/2), becomes small for small y. The very point is that it is not possible with 
totally inclusive methods to single out the photon-gluon fusion contribution. This would 
be possible by triggering on charm or bottom quarks, but those particles can only be 
detected indirectly, which is connected to losses in statistics. So, by combining both 
methods one can get a quite reliable method of extracting AG from DIS data. For 
the detection of the jets we use an algorithm which may be called a combined z-s-cone 
scheme [Eq. ([}])]. This is done in order to be able to work in the same jet scheme in 
Pepsi/Lepto and Mepjet because in this point the two codes differ substantially. The 
following parameters are set: 

0.04 

100 GeV 2 Pepsi 

100 GeV 2 Pepsi, Mepjet (2) 
1.0 

5 GeV . 

z min is the minimum of z = (P ■ Pj et )/(P ■ q) for the two jet four-momenta pj e t in a 
dijet event. The value z min = 0.04 is the default value of Lepto-6.5. s m j n is defined 
via s := (p + q) 2 , where p is the four-momentum of the incoming quark. The s cut is 
ineffective due to the large px cut. We furthermore require for the invariant mass of the 
two jets Sij := {pi +P2) 2 > 100 (GeV) 2 , where p\ and P2 are the two jet four-momenta. It 
has to be noticed that the cuts in Q 2 and pt are the smallest values possible in order to 
get reliable results within the LO cross sections. For smaller values NLO effects become 
dominant and the LO simulations are no longer applicable in a reliable way. -R m i n defines 



^min 
■^min 
"5ij',min 

PTmin 
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the minimal cone distance which two partons must have in order to belong to different 
jets. The cone distance is given by the Pythagorean sum of azimuthal deviation A0 and 
pseudo rapidity difference At/, i.e.: 



R = J(A v y + (A0) 2 . 



(3) 



Finally, for a^we use the one- loop expression with Aqcd — 151 MeV (which is the value 
for Aqcd over the bottom threshold), i.e.: 



a s (Q 2 ) 



Air 



11 - ln f ) ln(Q7A QCD ) 



(4) 



For a em we use the standard form as implemented in Jetset ||. For the interesting 
asymmetries, i.e., the ratio of the polarized and the unpolarized cross section, these 
coupling constants will cancel out. Regarding the parton distributions, we use the set 
'Gluon A' from Gehrmann and Stirling Jl0[ as the polarized one and the 'GRV leading 
order (LO) set' given in |~2] as the unpolarized one. For the observables discussed here 
we have checked that the qualitative behavior is the same for the standard LO scenario 



discussed in [FT]. The influence of several different models of AG on the polarized dijet 
cross section has been discussed in Iffl. 



We first focus on the angle between the spatial momenta of the exchanged pho- 
ton and the outgoing quark jet. From a theoretical point of view we can define an angle 
th as the angle between the exchanged photon and the outgoing quark jet in the cm 
system of the two outgoing jets. Then, in Fig. |2| the upward going jet 1 would be a 
quark jet and the downward going jet 2 would be an anti-quark jet or a gluon jet. Such 
a definition, although interesting for theoretical considerations, is not accessible experi- 
mentally because one cannot distinguish between quark, anti-quark and gluon jets. An 
experimentally accessible definition (0 exp ) would be that exp is the intersection angle 
(smaller than 90 degrees) between the line formed by the two jets in their cm system 
and the incoming photon. G th is connected with the variable z from Eq. (g) via the 
relation z = cos 2 (G th /2). In order to isolate the pure spin- dependent contribution one 
first regards the cross sections for anti-parallel (cfer(H)) and parallel spins (o?cj(t|)) of the 
colliding proton and electron beams. One then takes the mean and half of the difference 
of the differential cross sections of the two polarization configurations: 

da = l(da(U) + MW) > 

dAa = ~(da(U) -<fr(TT)) • (5) 

In Fig. | we compare the polarized and unpolarized th distributions and the predicted 
asymmetries as calculated with Pepsi (filled squares) and Mepjet (histogram). The 
large deviations for the photon-gluon fusion cross sections are due to the treatment of the 
quark masses and target remnant masses. For the shown Mepjet results quark masses 
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and target remnant masses are neglected completely. In Pepsi/Lepto they are taken 
into account via threshold effects only, i.e., the energy and momentum must be fulfilled 
when the quark masses and remnant masses are added to the jet four-vectors, whereas 
the cross sections in Pepsi/Lepto are calculated for massless partons. As Fig. |3| shows, 
both effects are nearly absent for the gluon-bremsstrahlung process, while they play a 
crucial role for the photon-gluon fusion. Moreover, the effect due to the target remnant 
mass is dominating. Effectively this remnant mass is included into the Lepto code in 
such a way that not only the massive quark anti-quark pair, but also the remnant has 
to be generated from the momentum four-vectors of the incoming photon and proton. 
In the case of the gluon bremsstrahlung such an effect is absent because as gluons are 
massless there is effectively no energy necessary to radiate a gluon. 
To demonstrate this we take an upgraded version of Mepjet, which takes into account 
the exact cross section for photon gluon fusion with quark masses and compare this 
to a Pepsi run where we substitute the target remnant by a simple massless gluon 
(see Tab. 1). In all cases we take the heavy quark masses to be m c = 1.35 GeV and 
mi, = 5 GeV, as used in standard Jetset. The two cross sections are qualitatively in 
agreement. Note that they may not agree exactly because Pepsi/Lepto takes quark 
masses into account only via threshold constraints, i.e., energy- momentum constraints 
inserted into massless cross sections, while Mepjet in the unpolarized case takes quark 
masses exactly into account. Unfortunately, there is so far no polarized version of the 
massive Mepjet available. Besides, as the agreement shows, the Pepsi/Lepto imple- 
mentation of quark masses is a quantitatively good approximation to the exact quark 
mass treatment, and, moreover, the absolute effect of the quark masses in the final state 
is only about 5-6% of the total cross section. When taking into account also the mass 
of the remnant, the photon-gluon fusion cross section is reduced by another 10%. There 
is no corresponding reduction in the bremsstrahlung cross section. This is due to the 
fact that for emitting a gluon, which is massless by itself, the energy can be in principle 
arbitrarily small. 

Moreover, the quark mass effects, as implemented in Pepsi enter the polarized and 
the unpolarized cross sections in the same way. For the asymmetries, displayed in the 
last row of Fig. ^|, the agreement between the two programs is consequently nearly exact. 
This on the other hand means that for asymmetry calculations we can use the massless 
program to get proper results. The numbers for the total cross sections are given in 
Tab. 1 for the unpolarized case and in Tab. 2 for the polarized case. The errors are 
the statistical MC errors, an inclusion of the systematic errors from parton distribution 
uncertainties etc. cannot be done in a reliable way. But taking the precision of present 
data as used in |IJ for example, one could accept a systematic error of at least 10% to 
be realistic. 
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cross sect, [pb] 
GB event 


cross sect, [pb] 
PGF-event 


MEPJET(l) massless version 


204.1 ±0.8 


1057.6 ±3.3 


Mepjet(2) exact quark masses 


205.8 ±1.1 


991.5 ±2.5 


Pepsi/Lepto(1) approx. quark masses 


206.4 ± 3.0 


980.4 ± 13.8 


Pepsi/Lepto(2) approx. quark 
and remnant masses 


205.5 ±3.0 


897.0 ± 12.7 



Table 1: Treatment of quark masses in Pepsi/Lepto and Mepjet in the unpolarized 
case. 

3 Asymmetry distributions and optimization of kine- 
matic cuts 

Fig. |] shows that the asymmetry is largest for angles th , where the axis formed by the 
two jets lies close to the spatial momentum of the incoming photon, and that the quark 
jets in the gluon bremsstrahlung are accumulated slightly for the unpolarized case and 
distinctly for the polarized case close to the incoming photon. This has quite a natural 
reason as the quark jet in gluon bremsstrahlung is ordinarily more energetic than the 
respective gluon jet. 

In Tab. 2 we give the integrated polarized cross sections. One should note again 
that the differences for the photon-gluon fusion come from the fact that in Mepjet 
quark masses are neglected. For a detailed discussion of the influence of quark masses 
on two and three jet rates see M and references therein. The most important question 
is how the significance of the asymmetry A = can be amplified by suitable cuts. We 
will discuss several kinematic variables in the following. Let us first concentrate on the 
angle O exp . As there is no possibility to distinguish with 100% reliability quark- and 
gluon jets, only the angle between the axis formed by the dijet in its cm frame and the 
incoming photon in the same frame exp is experimentally accessible. Of course then 
the asymmetry is large for small O exp , which can be already deduced from Fig. [| but on 
the other hand a cut 0^ p x also means losses in statistics. To get some feeling whether 
the signal could be improved, when imposing a O^ax cu ^> we define the significance sc 
by the negative logarithm to the basis ten of the relative error in the asymmetry A. 
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cross sect. \pb\ 
Pepsi(2) 


cross sect, [pb] 
Mepjet(I) 


A GB-event 
A PGF-event 


3.80 ± 0.12 
-37.71 ± 0.63 


3.93 ± 0.02 
-43.47 ± 0.15 



Table 2: Polarized cross sections (Aa) for the kinematic used in the paper, comparing 
Pepsi(2) and the Mepjet(I) program. 



(The choice of the basis ten is made to simplify rescaling for different luminosities.) 

S c:=-log 10 (<L4M) . (6) 

Moreover, the logarithmic representation has the advantage to make it easier to plot 
large numerical ranges for 5 A/ A. As the asymmetry is small for our kinematics, one can 
derive an approximative formula for the quadratic error of the asymmetry with only the 
luminosity and the unpolarized cross section: 



iV TT • N n 1 



5A = 2J ~ (7) 

V (*TT + ^Lfpb-Vunpol. 

We can eliminate the luminosity from the significance formula by defining the reduced 
significance scr: 



sc := —log 



10 



PAV (o- TT + a n f 



+ ilog 10 (L[pb- 1 ]) 

=: S cr + ilog 10 (P 2 L[pb- 1 ]) . (8) 

Here P = P e P p is the degree of polarization, which is the product of the degree of 
polarization of the proton beam P v and the electron beam P e , to include the effect of 
beam polarizations different from 100%. In this case the asymmetry is reduced while 
the absolute error, which depends approximately only on the unpolarized cross section, 
remains constant. In the following, if not something else is stated explicitely, the asym- 
metry plots will correspond to a degree of polarization of 100%. The luminosity L[pb _1 ] 
is given per relative polarization, i.e., = La^ and = La^. If sc is zero or 
even negative then the error is as large or even larger than the signal. For a realistic 
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experiment sc should at least yield a number of 0.6, then one can set four bins with a 
two sigma signal each. The definition of scr allows to find the optimal cut for a maximal 
significance by a look on a single diagram. It is just a very economical way of visualizing 
the optimal parameter without trying several MC runs. In the following we plot the 
asymmetries versus several variables together with the reduced significances scr for the 
corresponding minimal or maximal cuts. In order to demonstrate the usefulness of this 
concept we show in Fig. |] the asymmetry and scr for the experimentally accessible exp 
and z exp region. z exp is the minimum of the two z values that belong to the two jets in 
Fig. |2]. Both variables are closely related to each other because the more energetic one 
of the two parton jets is as compared to the other, i.e., the smaller z exp is, the closer 
this jet will be aligned to the incoming photon in the frame defined in Fig. ^|, i.e., the 
smaller is the angle G exp (z cxp = cos 2 ((180° — B exp )/2)). The asymmetry over both 
variables is displayed in Fig. ^. In both cases, the gluon bremsstrahlung reduces the 
asymmetry signal coming from photon-gluon fusion alone because the spin-dependent 
contributions to both processes differ in sign. From the shape of the asymmetries (left 
column in Fig. f|) one could suspect that a suitable Q^ p x or z^ p x cut would increase the 
measurable asymmetry signal. However, the significance plots on the right show that 
this is not the case. The reason for this is that the spin-dependent contribution Act is 
small, at most ten percent of the unpolarized cross section, and therefore the error [see 
Eq. (0)] is approximately anti-proportional to the square root of the unpolarized cross 
section alone. With decreasing 9^ x or z^ p x , as can be seen in Fig. ||, the unpolarized 
cross section is reduced and therefore the error increases, over-compensating the gain 
in the asymmetry signal. As an illustration Fig. [| shows the expected experimental 
errors for a luminosity of L = 100 pb _1 per relative polarization, i.e., L = 200 pb _1 in 
total. The assumed electron and proton polarization is 70% each. The asymmetry is 
sizable using realistic cuts for the HERA machine. The envisaged luminosity allows to 
extract a clear and distinct signal in five bins, for example. The first bin corresponds 
to O^jL = 28°. The reduced significance for this cut in Fig. |] is scr = —0.2232. For 
L = 100 pb _1 per relative polarization and a degree of polarization of P = 0.5 this 
corresponds to a relative error of 0.334. 

We now turn to the pt distributions in the lab system, see Fig. |6], where the maxi- 
mum and the minimum pt of the two jets is plotted. For pure photon-gluon fusion, 
the asymmetry is growing for larger pT,max- Taking into account gluon bremsstrahlung 
the asymmetry is not only reduced, but also decreases above pr,max = 25 GeV. The 
behavior for pT,min is similar. Moreover, as a glance to the reduced significances shows 
the relative error is increasing for increasing pt cuts so that we recommend to choose, 
in spite of the behavior of the asymmetries, the pt cut not larger than 5 (GeV) 2 . 

The same is true for the variable y, see Fig. |7j, where again by increasing the y m in 
cut one gains in asymmetry but loses significance. The only way to increase the signifi- 
cance is by choosing a definite minimal cut for the invariant dijet mass s^, which should 
be of the order s^mm = 250 GeV 2 . However, it was shown in [|J that for unpolarized 
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jet production it is advisable to choose Sj^min = 500 GeV 2 in order to minimize next to 
leading order effects. A corresponding analysis for the polarized case awaits the imple- 
mentation of the polarized NLO cross sections into Mepjet. As can be seen in Fig. |7|, 
the losses in statistics for such cuts are not considerable. 

The last variable discussed here is Q. The asymmetry distribution is nearly constant 
as a function of Q, so this is a completely uninteresting variable for improving the sig- 
nificance. On the other hand, as we show below, a suitable cut in Q allows to change 
the ratio of photon-gluon events to gluon bremsstrahlung events. In Fig. § we show the 
ratio of the total asymmetry and the pure gluonic asymmetry. The influence of gluon 
bremsstrahlung on the asymmetry increases for smaller Q^ in . The significance is con- 
siderably larger for small values of <3^ in , so one ma y think of choosing Q^ in as small as 
Qmm = 5 GeV 2 in the envisaged HERA experiments for extracting the polarized gluon 
density. 

In summary we have found that the two programs Pepsi and Mepjet (in the massless 
version) agree for the predicted asymmetries, as the mass effects enter the polarized and 
unpolarized cross sections in Pepsi in the same way and therefore do not contribute 
to the asymmetry. For a polarized HERA the invariant dijet mass Sy, m in plays a very 
special role, as it is the only variable where the significance of the dijet asymmetry 
signal can be slightly improved by increasing the minimal cut. For the invariant dijet 
mass one should choose an Sy jm i n cut between 250 and 500 GeV 2 . In all other cases 
studied here, it is advisable to choose small minimal cuts for exp , z,pT,y,Q 2 in order 
to maximize the significance. However, the final choice of those cuts will be dictated by 
future studies comparing LO and NLO polarized calculations. We have shown that the 
influence of quark masses on the gluon bremsstrahlung cross section is negligible, while it 
is substantial for the photon-gluon fusion. A very dominant effect for the photon-gluon 
cross section may also arise from the remnant mass, which makes corrections on the 
10% level. We have furthermore reproduced the observation that for a polarized HERA 
the asymmetry signal is sizable and highly significant for the envisaged luminosity of 
100 pb _1 per relative polarization and 70% beam polarization. 
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Jet 1 




Figure 2: Definition of the angle G: From a theoretical point of view an angle can 
be defined by the angle between the spatial momenta of the exchanged photon and the 
outgoing quark in the cm system of the two outgoing jets. Then the upward going jet 1 is 
a quark jet and the downward going an anti-quark or a gluon jet. But such a definition 
is not accessible experimentally because one cannot distinguish between quark, anti- 
quark and gluon jets. An experimental accessible definition (0 cxp ) is that O cxp is the 
intersection angle (smaller than 90 degrees) between the line formed by the two jets in 
their cm system and the incoming photon line. 
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Ada/da, GB-event, U Ada/da, PGF-event, 



Figure 3: Influence of quark-mass effects to the unpolarized and polarized dijet cross sec- 
tions. We show the results of the massless Mepjet program (histogram) and the Pepsi 
code (full squares), where quark- mass effects and remnant-mass effects are included. 
The units are in pb/ degree. 
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Figure 4: Asymmetry for proton beam A v and reduced significance (scr) for the exp - 
and z exp distributions. The histogram on the left hand side shows the asymmetry for 
photon-gluon fusion only. The triangles show the real asymmetry, which contains gluon 
bremsstrahlung as well as photon-gluon fusion. The pictures on the right show the 
significance (scr) plotted against the 6^ x cut and z^ x cut, respectively. The data are 
simulated with Pepsi. 
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Figure 5: Experimentally accessible asymmetry versus exp . The data are taken from 
the Pepsi run. The full circles denote the total asymmetry while the empty circles 
are the asymmetry if the gluon bremsstrahlung process is switched off. The horizontal 
bars denote the binning width. In the plot we assumed that proton and electron beam 
polarizations are 70% each. The error bars show the expected experimental error for 
100 pb~ x per relative polarization. 
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Figure 6: Asymmetry A p and reduced significance (scr) for the pT,max- and pT,min dis- 
tributions. The histogram on the left hand side shows the asymmetry for photon-gluon 
fusion only. The triangles show the real asymmetry, where gluon bremsstrahlung to- 
gether with the photon-gluon fusion is taken into account. The pictures on the right 
show the significance versus PT,max and PT,min lower boundaries. The data are simulated 
with Pepsi. 
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Figure 7: Asymmetry A p and reduced significance (scr) for the y- and s^- distributions. 
Histogram: no gluon bremsstrahlung; triangles: full asymmetry. Data from Pepsi. The 
plots on the right show the significance as a function of y(min) and ^(min), respectively. 
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Figure 8: Ratio of the total asymmetry to the pure gluon initiated asymmetry (top) and 
reduced significance scr as a functions of <5(min) (bottom); PEPSI data. 
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